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Abstract. The present study takes the CoRoT-7b exoplanet as an analogue for massive terres- 
trial planets to investigate conditions, under which intrinsic magnetic fields could be sustained 
in liquid cores. We examine the effect of depth-dependent transport parameters (e.g., activation 
volume of mantle rock) on a planet's thermal structure and the related heat flux across the core 
mantle boundary. For terrestrial planets more massive than the Earth, our calculations suggest 
that a substantial part of the lowermost mantle is in a sluggish convective regime, primarily 
due to pressure effects on viscosity. Hence, we find substantially higher core temperatures than 
previously reported from parameterized convection models. We also discuss the effect of melting 
point depression in the presence of impurities (e.g., sulfur) in iron-rich cores and compare cor- 
responding melting relations to the calculated thermal structure. Since impurity effects become 
less important at the elevated pressure and temperature conditions prevalent in the deep interior 
of CoRoT-7b, iron-rich cores are likely solid, implying that a self-sustained magnetic field would 
be absent. 
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1. Introduction 

CoRoT-7b is probably the most prominent discovery of the CoRoT space mission as it 
is the first extrasolar planet below ten Earth masses with firm observational constraints 
on planetary mass and radius. Although the planet's exact mass is still subject to debate 
(e.g, Pont et al. 2010 Hatzes et al. 2010), a total mass of about five times that of the 



Earth seems to be most favored (Bruntt et al. 2010 Queloz et al. 2009). Combined 



with the measured radius of (1.58 ± 0.10) (Bruntt et al. 2010), this suggests an 
average compressed density of (7.2 ± 1.8) Mg m~'^ for CoRoT-7b, being consistent with 
a terrestrial bulk composition (Fig. [T]). In the following, we take CoRoT-7b as a type 
example for a terrestrial exoplanet to investigate the physical state and thermal structure 
of planetary interiors. The purpose of this study is to understand the necessary conditions, 
under which an intrinsic magnetic field could be sustained in liquid cores of massive 
terrestrial planets. 



2. Method 

We use a one-dimensional, four-layer structural model that is similar to the model 
approach disclosed on the occasion of the XXVIIth lAU General Assembly held in Rio 



de Janeiro, Brazil (Wagner et al. 2009). The interior structure of CoRoT-7b is obtained 



by solving the mass and energy balance equations in conjunction with an equation of state 



(EoS) for the internal density distribution. The Birch-Murnaghan EoS (Birch 1947) and 



the Rydberg-Vinet EoS (Vinet et al. 20101 are often criticized because of inconsistencies 



in respect to the physics of the thermodynamic limit (e.g., Stacey & Davis 2004 1. Hence 
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Figure 1. Mass-radius relations for 
differentiated exoplanets (dashed 
lines) and homogeneous, self-com- 
pressible spheres of water, silicate, 
and iron (solid lines) ranging from 
1 to 10 M^. The relative position 
of the Earth is indicated by the 
triangle. Whereas the gray ellipse 
denotes the previously reported 



range for to tal mass Mp (Queloz 



let al. 2009|) and plane tary radius 
Rp ( |Leger et al. 2009] ), the black 
ellipse represents the recently revised 
range for Rp = (1.58 ± 0.10) Re and 



20101 



(5.2 ± 0.8) Me ( [Bruntt et ah] 



we improve previous models (e.g., Valencia et al. 2007) by implementing a generalized 
Rydberg EoS (Stacey 2005), facilitating extrapolation to exceptionally high pressures 



( .Wagner et al. 2010[ ) 



We adopt a mixing length formulation (Sasaki & Nakazawa 1986) to self-consistently 



calculate the thermal state of planetary mantles. This is an improvement to our previ- 
ously proposed model, in which we have used this approach only for the upper mantle to 
simulate a lithosphere within the uppermost part (Wagner et al. 2009). The basic idea 



behind the mixing length concept is that internally generated heat is primarily trans- 
ferred by vertical motion of fluid parcels which will entirely loose their individuality after 



migrating across size-dependent characteristic length scales. Yamagishi fc Yanagisawa 
|2001| demonstrated the feasibility of that concept for modeling of heat transfer within 
terrestrial planet interiors. The viscosity determining the temperature profile is mod- 
eled in the framework of a temperature- and pressure-dependent Arrhenius viscosity law, 
considering diffusion and dislocation creep as the most important creep mechanisms in 



sihcate aggregates (Ranalli 2001). The pressure-induced reduction of activation volume 



with depth is approximated as a vacancy in the material (O'Connell 19771. Solidification 



temperatures for the silicate-dominated mantle are obtained from ab initio molecular dy- 
namics simulations (Belonoshko et al. 2005). To construct melting relations under core 



conditions for pure iron and iron-rich alloys containing impurities, we use a parameter- 



ization based on the theory of tricritical phenomena (Aitta 2010) up to 800 GPa and 



extrapolate according to the well-known Lindemann law (Ross 1969). 



3. Results and Discussion 

First, we present modeling results as obtained by using CoRoT-7b as case study and 
calculating self-consistently the thermal structure within planetary mantles and cores. 
Figure [2] (left-hand side) shows the radial temperature profile of two plausible CoRoT-7b 
models of 5 that differ in the iron core mass fraction assumed, representing Earth- 
like (i.e., 32.6 wt.-%) and iron-depleted (i.e., 15.0 wt.-%) bulk compositions. From a 
fixed surface value of 1270 K (Leger et al. 2009), temperature increases rapidly across 



a thin thermal lithosphere, where heat is transferred by conduction and convection is 
absent. A mostly adiabatic temperature rise is observed across the converting mantle 
and underlying core, resulting for the Earth-like and iron-depleted case, respectively, 
in temperatures of 5210 K and 5320 K at the core-mantle boundary (CMB) and up 
to 7560 K and 6710 K in the planet's center. Compared to parameterized convection 
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Figure 2. Interior structure of CoRoT-7b corresponding to different iron core mass fractions: 
(Left) Radial distribut ion of temperature. An Earth-sized model and a reference geotherm 
(Stacey & Davis 20081 are shown for comparison. (Right) Radial distribution of hydrostatic 



pressure. 



models (e.g., Valencia et al. 2006), we find generally higher CMB temperatures. This is 
caused by a steeper temperature gradient prevalent in the deep interior that is mainly 
attributed to the pressure-induced increase of viscosity at elevated lower mantle pressures 
and temperatures. To illustrate environmental conditions in the deep interiors of massive 
exoplanets, the corresponding pressure profiles for the two CoRoT-7b cases are shown in 
Fig.[2]right-sided. Whereas CMB pressures of CoRoT-7b are twice as large as the central 
pressure of the Earth, the central pressure is as much as five times the pressure at the 
center of the Earth, depending on the planet's iron core mass fraction. 

Next, we investigate conditions for the possible existence of liquid cores, which are pre- 
requisite for magnetic field generation. In Figure [s] pressure-temperature (P-T) profiles 
for two Earth-like CoRoT-7b models are shown together with the corresponding silicate 
and iron melting relations. It is seen that iron melting temperatures are generally lower 
than solidus temperatures of MgSiOa (post-)perovskite under CMB conditions. The tem- 
perature profile shown in black corresponds to the Earth-like model case illustrated in 
Fig. [2] and yields an activation volume of 1.1 cm^ mol""'^ at the CMB. Lower mantle 
and core materials are expected to be solid because of extremely high compression at 
relatively low temperature. A temperature increase of about 3000 K at the CMB would 
be required to facilitate melting of the outer core. The temperature profile shown in gray 
corresponds to the latter case, for which a larger activation volume of 1.8 cm'^ mol^^ at 
the CMB would be required. Hence, the lower mantle becomes stiffer with increasing ac- 
tivation volume and, as a consequence, local temperature will rise until intersection with 
the iron melting curve. Also notable in that respect is that the melting intervall in the 
presence of impurities (e.g., siderophile elements like sulfur) would vanish at extremely 
high pressures (blue curves), according to (Aitta 2010). Due to large pressures within 



the deep interior of a 5 M0 exoplanet such as CoRoT-7b, the addition of impurities will 
not modify much the P-T relation compared to that of a pure iron core. This suggests 
that impurities should only have a minor effect on core melting within massive terrestrial 



exoplanets (Wagner et al. 20101. 



4. Conclusions 

It is concluded that the physical state of the deep interior of massive terrestrial ex- 
oplanets is strongly dependent on mantle rheology. Impurities have only a minor effect 
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Figure 3. Influence of activation volume V* and core impurities (e.g., sulfur) on the physical 
state of matter within a 5 exoplanet. 

on core melting, which can be explained conceptually by the large compression of mat- 
ter under core conditions. Due to the large effect of pressure on melting, a pure iron 
core is expected to be solid and, therefore, a self-generated magnetic field should be 
absent on massive terrestrial exoplancts like CoRoT-7b. Nevertheless, liquid cores can- 
not completely be ruled out, but in that case substantially larger activation volumes 
(> 1.8 cm'^ mol~^) compared to values predicted by the vacancy approach would be 
required to initiate core melting. Furthermore, a pressure-induced sluggish convection 
is prevalent in the lowermost mantle and should influence mantle convection pattern as 
well as the thermal evolution of massive exoplanets. 
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